Spin-spin interaction and spin-squeezing in an optical lattice 
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We show that by displacing two optical lattices with respect to each other, we may produce 
interactions similar to the ones describing ferro-magnetism in condensed matter physics. We also 
show that particularly simple choices of the interaction lead to spin-squeezing, which may be used 
to improve the sensitivity of atomic clocks. Spin-squeezing is generated even with partially, and 
randomly, filled lattices, and our proposal may be implemented with current technology. 
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Simulation of quantum many-body problems on a clas- 
sical computer is difficult because the size of the Hilbert 
space grows exponentially with the number of particles. 
As suggested by Feynman the growth in computa- 
tional requirements is only linear on a quantum computer 
1^, which is itself a quantum many-body system, and 
such a computer containing only a few tens of quantum 
bits may outperform a classical computer. A quantum 
computer aimed at the solution of a quantum problem is 
expected to be easier to realize in practice than a general 
purpose quantum computer, because the desired solu- 
tion is governed by physical interactions which are con- 
strained, e.g., by locality In essence, such a quan- 
tum computer is a quantum simulator with the attractive 
feature that the experimentalist can control and observe 
the dynamics more precisely than in the physical system 
of interest. In this Letter we describe how atoms in an op- 
tical lattice may be manipulated to simulate spin-spin in- 
teractions which are used to describe ferro-magnetism in 
condensed matter physics. We also show that with a spe- 
cific choice of interaction we may generate spin squeezed 
states which may be used to enhance spectroscopic 
resolution e.g., in atomic clocks. 

In Refs. [HQl two different methods to perform a co- 
herent evolution of the joint state of pairs of atoms in 
an optical lattice were proposed. Both methods involve 
displacement of two identical optical lattices with respect 
to each other. Each lattice traps one of the two internal 
states |0) and |1) of the atoms. Initially, the atoms are in 
the same internal state |0), the two lattices are on top of 
each other and the atoms are assumed to be cooled to the 
vibrational ground state in the lattice. Using a resonant 
pulse the atoms may be prepared in any superposition 
of the two internal states. The lattice containing the |1) 
component of the wavefunction is now displaced so that 
if an atom (at the lattice site k) is in |1), it is transferred 
to the vicinity of the neighbouring atom (at the lattice 
site fc -|- 1) if this is in |0), causing an interaction between 
the two atoms. See Fig. 0. The procedures described 
in this Letter follow the proposal in Ref. Q where, the 
atoms interact through controlled collisions. Also the op- 
tically induced dipole-dipole interactions proposed in |^] 



may be adjusted to fit into this framework. After the 
interaction, the lattices are returned to their initial posi- 
tion and the internal states of each atom may again be 
subject to single particle unitary evolution. The total 
effect of the displacement and the interaction with the 
neighbour can be tailored to yield a certain phaseshift (p 
on the |l)fe|0)fe-|-i component of the wavefunction, i.e.. 
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where \a)k (a = or 1) refers to the state of the atom 
at the fc'th lattice site. In it is suggested to build a 
general purpose quantum computer in an optical lattice. 
Such a general computer requires two- atom gates, which 
may be accomplished through the dynamics in (0) and 
single atom control, which is possible by directing a laser 
beam on each atom. We shall show that even without al- 
lowing access to the individual atoms, the lattice may be 
used to perform a highly non-trivial computational task: 
Simulation of a ferro-magnet. 

Our two level quantum systems conveniently de- 
scribe spin 1/2 particles with the two states |0)fc and 
|l)fe representing \jm)k = |l/2,-l/2)fe and |l/2,l/2)fe, 
where states \jm)k are eigenstates of the jz^fc-operator 
jz,k\j'm')k — 'ni\jm)k [h = 1). The phase-shifted com- 
ponent of the wavefunction in Eq. (^ may be iso- 
lated by applying the operator {j^^k + l/2){j^^k+i - 1/2), 
and the total evolution composed of the lattice transla- 
tions and the interaction induced phaseshift may be de- 
scribed by the unitary operator e^*^* with Hamiltonian 
H = xijz,k + l/2)(jz,fc+i - 1/2) and time t = (j)/x- In 
a filled lattice the evolution is described by the Hamil- 
tonian H = xLfcOz.fe + l/2)(i^,fc-Hi - 1/2), and if we 
are only interested in the bulk behaviour of the atoms 
we may apply periodic boundary conditions, so that the 
Hamiltonian reduces to 



<k,l> 



(2) 



where the sum is over nearest neighbours. By appropri- 
ately displacing the lattice we may extend the sum to 
nearest neighbours in two and three dimensions. Hzz 
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coincides with the celebrated Ising-model Hamiltonian 
introduced to describe ferro-magnetism. Hence, by 
elementary lattice displacements we perform a quantum 
simulation of a ferro-magnet. 

A more general Hamiltonian of the type 

= 'X-^^-kJzJ + VJx.kjx.l + ^jy.kjy,! (3) 

<k,l> 

may be engineered using multiple resonant pulses and 
displacements of the lattices: A resonant 7r/2-pulse act- 
ing simultaneously on all atoms rotates the jz-operators 
into jj^-operators, e^^y /"^ j z,ke~^^ /'^ = jx.k- Hence, 
by applying 7r/2-pulses, in conjunction with the displace- 
ment sequence, we turn H^z into Hxx^ the second term 
in Eq. (|^). Similarly we may produce Hyy, the third 
term in Eq. (|3|), and by adjusting the duration of the 
interaction with the neighbours we may adjust the coef- 
ficients X, 77 and A to any values. We cannot, however, 
produce Hf by simply applying Hzz for the desired time 
t, followed by H^x and Hyy^ because the different Hamil- 
tonians do not commute. Instead we apply a physical 
implementation of a well-known numerical scheme: The 
split operator technique. If we choose short time steps, 
i.e., small phaseshifts (p in Eq. (^, the error will only be 
of order 0^, and by repeated application of Hzz, Hxx and 
Hyy, we may stroboscopically approximate Hf. 

For a few atoms the system may be simulated numer- 
ically on a classical computer. In Fig. |2| we show the 
propagation of a spin wave in a one-dimensional string 
of 15 atoms which are initially in the |l/2,— 1/2) state. 
The central atom is flipped ai t = and a spin wave 
propagates to the left and right. The figure shows the 
evolution of < jz^k > for all atoms, obtained by repeat- 
edly applying the Hamiltonians Hzz, Hxx and Hyy with 
X — V — ^ and periodic boundary conditions. Small time 
steps dt = 0.1x~^ result in a stroboscopic approximation 
almost indistinguishable from the results of a direct nu- 
merical integration of Hf. 

A host of magnetic phenomena may be simulated on 
our optical lattice: Solitons, topological excitations, two 
magnon bound states, etc. By pumping a fraction of 
the atoms into the |l/2,l/2) state, we may also per- 
form micro-canonical ensemble calculations Q for non- 
vanishing temperature. Other procedures for introducing 
a non- vanishing temperature are described in Ref. . 

We now show how to generate spin squeezed states 
using the same techniques as discussed above. Signals 
obtained in spectroscopic investigations of a sample of 
two level atoms are expressed by the collective spin op- 
erators Ji — '^kji,k, and their quantum mechanical un- 
certainty limits the measurement accuracy, and e.g., the 
performance of atomic clocks. In standard spectroscopy 
with N uncorrelated atoms starting in the 1 1/2,— 1/2) 
state, the uncertainties AJ^, = \/ {Jx) ~ (Jx)^ and AJy 
are identical, and the standard quantum limit resulting 



from the uncertainty relation of angular momentum op- 
erators 



{AJx)^{AJy)^>\<Jz/2>f 
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predicts a spectroscopic sensitivity proportional to 
1/\/N. Polarization rotation spectroscopy and high pre- 
cision atomic fountain clocks are now limited by this sen- 
sitivity 0,|lll. In g it is suggested to produce spin 
squeezed states which redistribute the uncertainty un- 
evenly between components like Jx and Jy in (^, so that 
measurements, sensitive to the component with reduced 
uncertainty, become more precise. Spin squeezing result- 
ing from absorption of non-classical light has been sug- 
gested |l^ and demonstrated experimentally Ref. 
[Q presents an analysis of squeezing obtained from the 
non-linear couplings H = x^x and H = x{Jx ~ Jy)- 
For neutral atoms, such a coupling has been suggested in 
the spatial overlap of two components of a Bose-Einstein 
condensate j]!^ . Spin squeezing in an optical lattice has 
two main advantages compared to the condensates: The 
interaction can be turned on and off easily, and the local- 
ization at lattice sites increases the density and thus the 
interaction strength. The product of two collective spin 
operators involve terms jx,kjx,i for all atoms k and /, and 
this coupling may be produced by displacing the lattices 
several times so that the 1 1/2, 1/2) component of each 
atom visits every lattice site and interacts with all other 
atoms. In a large lattice such multiple displacements are 
not desirable. We shall show, however, that substantial 
spin-squeezing occurs through interaction with only a few 
nearby atoms, i.e., for Hamiltonians 



H 



,kjx^l 
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and 



H — ^^XkA{jx,kjx,l - jy,kjy,l), 
k,l 



(5) 



(6) 



where the coupling constants Xk,i between atoms k and I 
vanishes except for a small selection of displacements of 
the lattices. 

Expectation values of relevant angular momentum op- 
erators and the variance of the spin operator Jg — 
cos{9)Jx + sm{0) Jy may be calculated for an initially 
uncorrelated state with all atoms in 1 1/2,— 1/2), prop- 
agated by the simple coupling (||). If each atom visits 
one neighbour Xk,i = X^k+i.i, we get the time dependent 



variance of the spin component J_j^ 
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The mean spin vector is in the negative z direction and 
has the expectation value 
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< 'h >= -—cos'ixt). 



(8) 



For small values of xt, AJ_^/4 decreases linearly with 
Xt whereas | < > | decreases proportional to (x*)^, 
hence AJ_7r/4 falls below | < Jz/2 > |, and the spin is 
squeezed. 

In Fig. ^ we show numerical results for 15 atoms in 
a one-dimensional lattice with periodic boundary condi- 
tions. Fig. 1^ (a) shows the evolution of (AJg)^ when 
we apply the coupling (|^) and visit 1, 2, and 3 neigh- 
bours. The squeezing angle 9 = —tt/A is optimal for 
short times 1- For longer times the optimal an- 

gle deviates from — 7r/4, and we plot the variance (AJg)^ 
minimized with respect to the angle 9. We assume the 
same phaseshift for all collisions, i.e., all non- vanishing 
Xk,i are identical. 

For spectroscopic investigations not only the variance 
of a spin component is relevant. In |^ it is shown that 
if spectroscopy is performed with N particles, the reduc- 
tion in the frequency variance due to squeezing is given 
by the quantity 



NjAJg 



(9) 



In Fig. ^ (b) we show the minimum value of obtained 
with the couplings (|^) and (0) as functions of the number 
of neighbours visited. Fig. (b) shows that the coupling 
(3) produces better squeezing than (||). The couphng 
(5), however, is more attractive from an experimental 
viewpoint. Firstly, all jx.k operators commute and we do 
not have to apply several displacements with infinitesi- 
mal durations to produce the desired Hamiltonian. We 
may simply displace the atoms so that they interact with 
one neighbour to produce the desired phaseshift 0, and 
then go on to interact with another neighbour. Secondly, 
if the jx,kjx,i coupling involves a phaseshift 0, the oper- 
ator —jy.kjy,i requires the opposite phaseshift —ip. This 
requires a long interaction producing 27r — (/>, or a change 
of the interaction among the atoms, i.e., a change of the 
sign of the scattering length in the implementation of [p| . 

Like the analytic expression for x^ obtained from , 
the results shown in Fig. || (b) are independent of the to- 
tal number of atoms as long as this exceeds the "number 
of neighbours visited" . When all lattice sites are visited 
we approach the results obtained in Ref. i.e., a vari- 
ance scaling as N^^^ and a constant for the couplings (|^) 
and d). 

So far we have assumed that the lattice contains one 
atom at each lattice site and that all atoms are cooled 
to the vibrational ground state. The present experimen- 
tal status is that atoms can be cooled to the vibrational 
ground-state, but with a filling factor below unity |jl^. A 
mean filling factor of unity is reported in , but when 
at most a single atom is permitted at each lattice site 



a mean occupation of 0.44 is achieved. It has been sug- 
gested that a single atom per lattice site may be achieved 
by filling the lattice from a Bose-Einstein condensate [0 . 

To describe a partially filled lattice it is convenient 
to introduce stochastic variables hk, describing whether 
the fc'th lattice site is filled ~ 1 or empty hk — 0. 
The interaction may be described by the Hamiltonian 
H = J2k,i Xk,ihkijz,k + l/2)hi{jz^i - 1/2), where the sum 
is over all lattice sites k and I. If we, rather than just 
displacing the atoms in one direction, also displace the 
lattice in the opposite direction, so that Xk,i is symmetric 
in k and Z, we may produce the Hamiltonian 



H y^^Xk,ihkjx,khijx 



(10) 
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This Hamiltonian models ferro-magnetism in random 
structures, and it might shed light on morphology prop- 
erties, and, e.g., percolation [Q. Here we shall restrict 
our analysis to spin-squeezing aspects, since these are 
both of practical interest, and they represent an ideal 
experimental signature of the microscopic interaction. 

In Fig. H we show the result of a simulation of squeez- 
ing in a partially filled one dimensional lattice. Each 
lattice site contains an atom with a probability p, and 
the size of the lattice is adjusted so that it contains 15 
atoms. In Fig. ^ (a) we show the decrease in the variance 
of Jg, averaged over 20 realizations and minimized with 
respect to 9. Lines indicate the predictions from the time 
derivatives at i = 



^(AJ_^/4)2 = -l^xk,i < hkhi > 
d 



k,l 



dt 



Jz=0, 



(11) 



where < h^hi > denotes the ensemble average over the 
distribution of atoms in the lattice, i.e., the two atom 
correlation function. In Fig. ^ (b) we show the minimum 
value of for different filling factors p as a function of 
the number of neighbours visited. The calculations con- 
firm that even in dilute lattices, considerable squeezing 
may be achieved by visiting a few neighbours. 

In conclusion we have suggested a method to simu- 
late condensed matter physics in an optical lattice, and 
we have shown how the dynamics may be employed to 
produce spin-squeezing. We emphasize the moderate ex- 
perimental requirement for our scheme. With the two 
internal states represented as hyperfine structure states 
in alkaline atoms, all spin rotations may be performed 
by Raman or RF-pulses acting on all atoms simultane- 
ously, and lattice displacements may be performed by 
simply rotating the polarisation of the lasers Q. With 
the parameters in , the duration of the sequence in Fig. 
|l| can be as low as a few micro-seconds. Following our 
suggestion spin-squeezing may be produced in dilute op- 
tical lattices, and implementation is possible with current 
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technology. The resulting macroscopic decrease in pro- 
jection noise has several promising applications in tech- 
nology and quantum physics, and it provides an experi- 
mental signature of the microscopic interaction between 
the atoms. 
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FIG. 1. (a) Two overlapping lattices trapping the two inter- 
nal states |0) (black circle) and |1) (white circle). By resonant 
laser pulses the atoms can be prepared in any superposition 
of the two internal states, (b) The lattices are displaced so 
that if an atom is in the |1) state, it is moved close to the 
neighbouring atom if this is in |0) causing an interaction be- 
tween the two atoms, (c) The lattices are returned to their 
initial position and the atoms may be prepared in new super- 
positions by resonant pulses. 
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FIG. 2. Propagation of a spin wave in a one dimensional 
string. The central atom is flipped at t = 0, and repeated 
application of H^z, H^x and Hyy results in a wave propagat- 
ing to the left and right. The figure shows the evolution of 
< jz,k > for all atoms (fc). 
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FIG. 3. Squeezing in a one-dimensional lattice with 15 
atoms, (a) Evolution of (AJg)^ with the coupling (^ and 
interaction with 1, 2, and 3 neighbours (full, dashed, and 
short dashed line, respectively) . (b) The optimal value of the 
squeezing parameter obtained with the coupling (^) (o) 
and (^) (4-). Lines are shown to guide to the eye. 
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FIG. 4. Spin squeezing in a partially filled one dimensional 
lattice containing 15 atoms, (a) Evolution of (AJg)'^ in a lat- 
tice with a filling factor p = 50% and displacements to 1, 2, 
and 3 neighbouring sites (full, dashed, and short dashed curve 
respectively). Dotted lines represent the predictions from Eq. 
(pT). (b) Minimum attainable squeezing parameter for fill- 
ing factors p=100% (o), 50% {+), 25% (□), and 10% (x) as 
functions of the number of sites visited. 
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